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TE/TM  Height  Profile  Measurements  of 
Transpolar  VLF  Signals 


1.  INTRODU  CTIOIN 

Because  low  and  very  low  frequency  (LF/VLF)  radio  signals  propagate  with 
relatively  low  attenuation  within  the  Earth-ionosphere  waveguide,  they  are  used  ex¬ 
tensively  for  very  long  range,  survivable  communications.  Conventional  LF  /VLF 
systems  employ  large,  vertical  transmitting  antennas  that  radiate  transverse  mag¬ 
netic  (TM)  waves.  Recently,  it  has  become  practical  to  radiate  LF/VLF  energy 
from  the  aircraft,  the  TWAs  are  unspooled  to  lengths  approaching  half-wave  reso¬ 
nance,  making  them  very  efficient  radiators.  Due  to  air  drag,  the  TWAs  are 
nearly  horizontal,  so  they  couple  to  transverse  electric  (TE)  modes  as  well  as  to 
the  TM  modes.  TE  signals  have  fields  that  tend  to  be  very  weak  near  the  ground 
and  strong  at  aircraft  altitudes,  while  TM  signals  tend  to  be  strongest  near  the 
ground.  Thus,  there  is  increasing  interest  in  exploiting  polarization  diversity  (the 
combination  of  TE  and  TM  modes)  to  enhance  air-to-air  capabilities  of  the  Minimum 
ation  of  TE  and  TM  modes)  to  enhance  air-to-air  capabilities  of  the  Minimum 
Essential  Emergency  Communications  Network  (MEECN). 

Although  well -understood  theoretically,  TE  signals  have  not  been  measured 
nearly  as  extensively  as  those  of  the  more  commonly  used  TM  modes.  In  the  mid- 
1970s,  the  Rome  Air  Development  Center  (RADC)  conducted  a  series  of  balloon 
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probe  experiments  in  which  TE  and  TM  field  strength  profiles  of  a  number  of 

ground-based  and  airborne  LF/VLF  transmitters  were  measured  to  altitudes  of 

about  20km.  *  RADC  also  equipped  a  U-2  aircraft  to  measure  how  TE  and  TM  sig- 

2 

nals  varied  with  range  from  airborne  transmitters.  In  1979,  the  first  of  a  series 
of  Super  ARCAS  rocket  probes  carrying  LF/VLF  sensors  was  launched  from  the 
MASA  Wallops  Flight  Center,  Wallops  Island,  Virginia,  to  sample  the  TE  and  TM 

3 

noise  environment.  Subsequent  launches  observed  height  profiles  of  TE/TMsig- 

4 

nals  from  airborne  terminals  at  long  distances  over  sea  water.  A  comprehensive 
review  of  the  results  from  these  and  other  experiments  was  presented  at  a  meeting 
of  the  Electromagnetic  Wave  Propagation  Panel  of  AGARD  in  Brussels,  Belgium, 
in  1981. 5 

With  the  continuing  development  and  upgrade  of  low  frequency  components  of 
MEECN,  including  the  modernization  of  TWA  transmitters  on  the  Airborne  Com¬ 
mand  Posts  (ABNCP)  as  well  as  the  incorporation  of  Miniature  Receive  Terminals 
(MRT)  on-board  strategic  bombers,  an  adequate  database  pertaining  to  LF/VLF 
signal  propagation  and  noise  in  the  polar  regions  at  aircraft  altitudes  becomes 
critically  important.  Of  special  concern  are  the  characteristics  of  TE  and  TM 
modes  propagating  over  poorly  conducting  arctic  regions,  under  both  normal  and 
disturbed  ionospheric  conditions.  Under  sponsorship  of  the  Defense  Nuclear  Agency 
Agency  (DNA)  and  the  Electronic  Systems  Division  (ESD),  RADC  conducted  rocket 
experiments  in  1983,  using  standard  ARCAS  rockets  launched  from  Thule  AB, 
Greenland,  to  measure  height  profiles  of  TE  and  TM  signals  and  noise  within  the 
polar  cap  to  altitudes  of  about  60  km. 

This  report  presents  the  results  obtained  from  a  rocket  flight  launched  at 
0600  UT  on  21  September  1983  under  quiet  ionospheric  conditions.  The  payload 
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and  ground-based  instrumentation  are  described,  the  trajectory  and  attitude  orien¬ 
tation  compensations  applied  to  the  raw  data  are  derived,  and  TM  and  TE  height 
profiles  for  a  number  of  ground  and  airborne  LF/VLF  transmitters  are  presented 
and  discussed. 


2.  INSTRUMENTATION 

The  rocket  payload  (Figure  1)  was  designed  with  two  separate  receiving  chan¬ 
nels,  to  measure  orthogonal  magnetic  field  components  of  incident  waves.  Space 
considerations  required  the  antenna  mounted  transverse  to  the  rocket  axis  (the 
"TM  channel")  to  be  constructed  of  two  separate  10.2  cm-long  ferrite  cores,  each 
with  350  turns,  connected  in  series,  while  the  longitudinal  antenna  (the  "TE  chan¬ 
nel")  was  made  of  one  continuous  20.  3  cm  length  of  ferrite  with  650  turns.  A  fiber¬ 
glass  nosecone  covered  the  antennas  during  flight.  The  difference  in  physical 
length,  and,  therefore,  the  sensitivity,  between  the  two  antennas  was  calibrated 
by  placing  the  payload  (without  the  ARCAS  body)  in  a  radio  frequency  magnetic 
field  of  known  strength. 

Each  antenna  was  connected  to  a  high-gain,  broadband  (15  to  50kHz,  3  dB 
bandwidth)  amplifier  whose  output  frequency -modulated  a  separate  sub-carrier 
oscillator  (VCO  in  Figure  2).  The  two  sub-carrier  signals  then  were  mixed  to 
form  a  composite  that,  in  turn,  frequency-modulated  a  lOOmW  S-Band  (2252.5 
MHz)  telemetry  transmitter  with  its  associated  slot  antenna.  A  1  MHz  oscillator 
signal  also  was  mixed  into  the  composite  signal,  to  provide  a  method  for  estimat¬ 
ing  slant  range  to  the  rocket  payload  during  its  flight. 

The  telemetry  receiving  system  and  associated  ground-based  instrumentation 
are  shown  in  Figure  3.  At  Thule  AB,  the  telemetry  antenna,  S-band  receiver, 

FM  discriminator,  and  tape  recorder  were  supplied  by  Det  3  of  the  Air  Force  Sat¬ 
ellite  Control  Facility  (AFSCF),  The  video  output  of  the  telemetry  receiver  cor¬ 
responded  to  the  composite  signal  (sub-carriers  plus  1  MHz  sine  wave)  present  at 
the  input  of  the  rocket  transmitter.  Replicas  of  the  original  TE  and  TM  signals 
present  at  the  output  of  the  broadband  amplifier,  as  well  as  the  1  MHz  oscillator 
signal  in  the  rocket,  were  obtained  by  filtering  and  processing  through  FM  dis¬ 
criminators;  these  replicas  were  recorded  on  separate  analog  tape  tracks  to  pre¬ 
serve  the  bandwidth  of  the  two  rocket  data  channels  and  allow  post -flight  proces¬ 
sing  in  different  ways,  depending  on  the  study  requirements. 

To  estimate  the  slant  range  and  thus  derive  altitude,  the  1  MHz  replica  was 
input  to  the  "A"  channel  of  a  two-channel  oscilloscope,  while  another  1  MHz  signal 
derived  from  a  local  source  was  applied  to  the  "B"  channel.  In  the  "A  +  B"  mode, 
the  oscilloscope  sums  the  two  signals;  the  two  oscillators  are  brought  into  syn- 
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Figure  2.  Rocket  Instrumentation  for  Measuring  TM  and  TE  Height-Gain  Profiles. 
The  1  MHz  oscillator  provides  a  signal  for  estimating  slant  range  from  the  tele¬ 
metry/tracking  antenna 


chronism  by  adjusting  the  phase  of  the  ground-based  source  to  give  a  steady  pre¬ 
sentation  while  the  rocket  is  still  on  the  ground  in  the  launcher.  During  flight,  the 
phase  of  the  1  MHz  replica  lags  that  of  the  local  oscillator  due  to  the  increase  in 
distance;  by  counting  the  cycles  of  phase  shift,  the  change  in  distance  could  be  de¬ 
termined. 

A  tunable  narrowband  receiver  with  a  vertical  loop  antenna  was  also  operated 
in  conjunction  with  the  rocket  flight.  This  was  used  to  measure,  at  ground  level, 
the  signals  from  several  of  the  transmitters  monitored  by  the  rocket -borne  re¬ 
ceivers. 


3.  TRAJECTORY  AND  ATTITUDE  ANALYSIS 
The  ARCAS  flight  trajectory 

T(x,  y,  z)  =  fx(t),  y(t).  z(t)]  (1) 

can  be  determined  from  measurements  of  azimuth,  cp  ,  and  elevation,  0  ,  supplied 
by  AFSCF  in  printout  form,  combined  with  estimates  of  slant  range,  R,  obtained 
from  a  count  of  cycle  slips  between  the  1  MHz  oscillator  on  the  rocket  and  a  local 
1  MHz  oscillator.  For  1  MHz,  a  change  of  one  cycle  corresponds  to  a  change  in 
slant  range  of  0.3  km.  The  Cartesian  coordinates  can  be  found  using 
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Figure  3.  Ground-Based  Instrumentation  for 
Recording  TM  and  TE  Height -Gain  Profiles. 
The  rocket  trajectory  is  computed  using  azi¬ 
muth  and  elevation  from  the  tracking  antenna 
in  combination  with  range  estimated  from  the 
rocket  1  MHz  phase  lag 


x(t)  =  R  cos  0  cos  cp 
y(t)  =  R  cos  0  sin  cp 

z(t)  =  R  sin  0  (2) 

where  the  origin  of  coordinates  is  the  tracking  station  (see  Figure  4). 

A  sounding  rocket  trajectory  approximates  a  parabola  in  the  plane  determined 
by  the  burn-out  coordinates,  the  apogee,  and  the  impact  point.  In  the  x-z  olane. 
a  parabolic  trajectory  may  be  expressed 

z  -  b  =  -  (l/2p)  (x  -  a)2  (3) 

where  [a,  b]  are  the  apogee  coordinates  and  p  is  a  constant  known  as  the  semilatus 
rectum.  There  are  about  100  measured  trajectory  data  points  [x^  =  x(t.),  =  z(t.)] 
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Figure  4.  Coordinate  System  for  the  Rocket  Trajectory 

between  burn-out  and  apogee.  This  data  is  used  to  find  a  "best  fit"  parabola  by 
calculating 

A.  =  <l/2p)  =  -<z.  -  b')/(x.  -  a')2  (4) 

where  [a1,  b ' ]  is  an  estimate  of  the  apogee  coordinates,  then  computing  the  average 

<A>  =  (  IA.)/N,  (5) 

and  finally,  forming  the  mean  square  difference 

MSD  =  [  I(A.  -  <A>)2]/N  (6) 

This  process  is  repeated,  using  new  estimates  for  [a1,  b'],  until  coordinates  [a,  b] 
are  found,  which  makes  the  MSD  a  minimum.  Then  <A  >  and  [a,  b]  are  the  "best 
fit"  constants  for  the  parabola: 

z  =  -  <A>  (x  -  a)^  +  b,  (7) 

There  was  no  explicit  measurement  of  the  rocket  orientation,  although  it  is 
crucial  to  determining  the  fraction  of  the  TE  and  TM  fields  received  on  each  ortho¬ 
gonal  loop.  Instead,  it  was  assumed,  since  the  rocket  did  not  leave  the  atmos- 
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phere,  that  the  fins  would  constrain  the  rocket  axis  to  coincide  with  the  total  veloc¬ 
ity  vector,  V,  tangent  to  the  trajectory.  The  flight  elevation  angle  can  be  calcu¬ 
lated  from  Eq.  (7): 

Z  =  tan-1  [-  2<A>  (x  -  a)]  (8) 

where  Z  is  the  angle  between  the  horizontal  and  the  total  velocity  vector. 

Eqs.  (7)  and  (8)  express  rocket  altitude  and  flight  elevation  angle  as  functions 
of  x,  the  horizontal  distance.  However,  the  trajectory  data  are  obtained  as  a 
function  of  time,  that  is,  R  =  R(t),  0  =  0  (t),  and  cp  =  cp(t).  The  relation  between 
x  and  t  for  a  parabola  is  simply 

x(t)  =  Bt  +  C  (9) 

The  constants  B,  C  can  be  determined  by  linear  regression  using  the  measured 
values  [x.,  t.].  Then 

l  l 

z(t)  =- <A>(Bt  +  C  -  a)2  +  b  (10) 


C(t)  =  tan"1  [-  2  <A> (Bt  +  C  -  a)] 


Eq.  (10)  can  be  used  to  transform  the  TE  and  TM  channel  output  voltages  to  func¬ 
tions  of  altitude. 

The  fraction  of  a  "pure"  TE  or  TM  mode  actually  received  on  each  channel 
can  be  determined  using  Eq.  (11)  in  conjunction  with  the  source  azimuth,  E,  ,  and 
the  azimuth  of  the  rocket  trajectory,  cp .  The  TM  magnetic  field  vector  from  a  TM 
source  at  azimuth  E,  is 

Hm  =  HQmfsin  (  (p  -  C  )  x  +  cos  (<p  -  C)  y]  (12) 

where  Hgm  is  the  magnitude  of  the  incident  field  (in  ampere/meter).  To  find  the 
fraction  of  the  total  TM  field  in  the  direction  of  the  TM  antenna  axis,  denoted  as 

A 

A  ,  we  construct  a  unit  vector  b  perpendicular  to  the  rocket  axis: 
mm  r 

A  A  A  As 

b  =  -  cos  co  t  sin  t;  x  +  sin  u>  t  y  +  cos  to  t  cos  Z  z,  (13) 

s  s  s 

where  cog  is  the  angular  spin  rate  of  the  rocket,  then  take  the  dot  product  of 


>  *.|l' 


i, 

i 


Eqs.  (12)  and  (13): 


(Hm  •  bl/H^  =  -  cos  o)gt  sine  sin(  cp  -  £  )  +  sino>gt  cos(  cp  -  C  ) 

The  spin  modulates  the  magnitude  of  this  component.  The  rms  value  is 


A _ =  [(ft  *b)  J/H, 

mm  m  rms  Om 


^1  -  sin2  ( cp  -  O  cos2  e 


A  unit  vector  in  the  direction  of  the  TE  antenna  axis  (along  the  rocket  axis)  is 
given  by 

£  =  cose  x  +  sine  z  (16) 

so  the  fraction  of  TM  field  in  direction  to  is  the  dot  product 

Ame:^m“)/HOm  =sin(<P'^cosC  (17) 

The  magnetic  field  from  a  TE  source  has  both  vertical  and  longitudinal  com¬ 
ponents,  but  everywhere  except  very  close  to  the  ground,  the  vertical  component 
dominates: 

He  =  HQe£  (18) 

where  Hq6  is  the  magnitude  of  the  incident  (vertical)  field.  The  fraction  in  the  TM 
antenna  axis  direction  is 


(H  •  b)/H  =  costo  t  cos  e 

6  00  S 


which  varies  due  to  spin. 

The  effective  (rms)  value  is 


-*•  ~  ,  cose 

A  =  [(H  •  b)  ]/Hn  =  — ~ 

em  1  e  rms'  Oe  ^ 


In  the  TE  antenna  axis  direction,  the  fraction  is 


Aee  =  (He  *  “)/H0e  =  sin  ? 


4.  DATA  REDUCTION 


At  the  conclusion  of  the  launch  campaign,  there  existed  an  analog  data  tape  con¬ 
taining  broadband  LF/VLF  signal  voltages  from  the  TM  and  TE  receivers  as  func¬ 
tions  of  time.  The  first  step  in  data  reduction  was  to  prepare  digital  tapes  for 
each  frequency  of  interest.  Under  the  control  of  a  Tektronix  4054  computer,  a 
Nicolet  Model  446B  Spectrum  Analyzer  performed  FFTs  on  the  broadband  analog 
voltages;  by  first  setting  the  analyzer  bandwidth  to  50  kHz,  numerous  signals  usu¬ 
ally  could  be  observed  on  the  400 -point  display,  in  most  cases  corresponding  to 
known  ground-based  or  airborne  sources.  Next,  using  higher  resolution,  a  2  kHz 
bandwidth  was  displayed,  centered  on  a  frequency  of  interest.  We  selected  20 
points  (bandwidth  =  100Hz)  for  the  noise  and  5  points  (25 Hz)  for  the  signal,  aver¬ 
aged  each  set  of  points,  and  recorded  the  averages  once  per  second  for  the  duration 
of  the  4-min  flight.  A  flow  chart  of  the  data  reduction  steps  is  shown  in  Figure  5. 

The  next  step  was  to  plot  and  compare  the  signal  and  noise  averages  to  deter - 
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Figure  5.  Data  Reduction  Flow  Chart  for  Mode  Profile  Analysis 


mine  whether  an  adequate  signal -to-noise  ratio  (S/N)  existed  for  the  measurements. 
We  assumed  that  if  S/N  >  6  dB,  we  could  process  the  signal  data  without  having  to 
correct  for  noise  contamination. 

The  signals  in  each  channel  actually  were  due  to  a  combination  of  TE  and  TM 
fields  because  of  rocket  orientation  and  spin.  Since  the  rocket  was  only  slightly 
tilted  from  vertical  for  most  of  the  upward  portion  of  the  flight,  we  assumed  as  a 
first  approximation  that  each  channel's  signal  was  unaffected  by  the  orthogonal 
polarization.  Initial  estimates  of  TE  and  TM  field  strengths  then  were  made  sim¬ 
ply  by  applying  Amm  and  Age  to  adjust  for  the  small  tilt.  To  test  the  validity  of 
these  estimates,  we  calculated  the  cross-coupling  expected  in  each  channel,  using 
the  initial  estimates  and  applying  the  appropriate  cross-geometric  factors  Amg 

and  A  .  We  assumed  that  the  first  approximation  was  indeed  a  good  estimate  of 
em 

the  field  so  long  as  the  calculated  cross -coupling  was  6  or  more  dB  smaller  than 
the  actual  recorded  channel  voltage.  Where  this  was  not  the  case,  the  first  ap¬ 
proximation  was  considered  to  be  invalid,  and  we  made  no  claim  for  a  measure¬ 
ment  of  that  component. 


5.  SIGNAL  SOURCES,  LOCATIONS,  AND  PROPAGATION  PATHS 

Within  the  15  to  50kHz  bandwidth  of  the  rocket  receivers,  14  VLF/LF  signals 
were  detected  (S/N  >  6dB)  on  at  least  one  channel  over  some  portion  of  the  flight. 
Thirteen  of  these,  found  between  15  and  27  kHz,  can  be  seen  in  Figure  6;  this  com¬ 
posite  spectrum  was  obtained  by  repeatedly  playing  the  analog  tape  and  selecting 
2  kHz  spectra  viewed  on  the  analyzer  when  the  signals  appeared  to  be  the  strongest. 
Due  to  the  spin  of  the  TM  antenna,  the  spectra  of  sufficiently  narrow  signals  (band¬ 
width  <  spin  frequency)  are  observed  to  be  split  into  two  lines,  each  with  half  the 
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Figure  6.  Composite  Spectrum  of  Received  TM  Signals 


energy  of  the  original  line.  On  the  other  hand,  signals  whose  bandwidths  are 
much  wider  than  the  spin  frequency,  that  is,  "spread  spectrum"  signals,  show 
only  a  small  increase  in  width  since  only  lines  at  or  near  the  edges  are  displaced 
outside  the  original  spectrum.  Of  the  14  detected  signals,  13  tentatively  were 
identified  by  frequency  as  specific  LF/VLF  transmitters;  their  propagation  paths 
to  Thule  are  shown  in  Figure  7,  and  certain  relevant  parameters  are  given  in 
Table  1. 

All  the  identified  sources  except  TACAMO  used  ground-based  vertical  anten¬ 
nas.  The  TACAMO  aircraft  was  flying  at  20,  000  ft  over  the  western  Atlantic,  on 
a  heading  approximately  perpendicular  to  the  azimuth  of  Thule  and  at  sufficiently 
high  speed  to  cause  its  trailing-wire  antenna  to  assume  a  nearly  horizontal  orien¬ 
tation  broadside  to  the  TACAMO-to-Thule  propagation  path. 


Figure  7.  Propagation  Paths  to  Thule.  On  21  September,  the  day/night  terminator 
passed  through  the  North  Pole  (NP) 


Table  1.  Transmitter  Listing  by  Frequency 


Measured 

F  requency 
(kHz) 

Source  Identification 

Call  Location 

A  zimuth 
at  Thule 

Propagation 

Distance 

(kM) 

ERP 

(kW) 

15.  9 

LY 

Bordeaux,  France 

098 

4640 

33 

16.4 

JXZ 

Aldra,  Norway 

068 

2810 

[NA] 

16.  8 

FUB 

St.  Assisse,  France 

093 

4260 

30 

17.  1 

UMS 

Gorki,  U.  S.  S.  R. 

048 

4530 

315 

17.  8 

NAA 

Cutler,  Maine 

178 

3550 

1000 

18.  15 

[Unknown] 

19.05 

GQD 

U.K. 

096 

3740 

[NA] 

19.6 

GBZ 

Criggion,  U.K. 

096 

3740 

10 

20.25 

IDO 

Rome,  Italy 

086 

5260 

[NA] 

21.4 

NSS 

Annapolis,  Md. 

190 

4180 

266 

23.4 

NPM 

Haiku,  Hawaii 

275 

7670 

530 

24.  8 

NLK 

Jim  Creek,  Wash. 

248 

3900 

234 

26.  1 

TACAMO 

38°N,  72 °W 

184 

4280 

100 

48.  5 

[NA] 

Silver  Creek,  Nebr. 

217 

4140 

50 

6.  HEIGHT  PROFILES 

In  Figures  8  to  16,  height -gain  profiles  of  9  of  the  14  detected  sources  are 
shown.  In  each  of  the  figures: 

(1)  The  abscissa  is  given  in  dB  relative  to  1  uV/m,  indicating  that  voltage-to- 
field  strength  calibrations  have  been  applied  (the  spin  of  the  TM  antenna  has  been 
taken  into  account  by  multiplying  the  TM  signal  strength  by  V2). 

(2)  In  (a)  and  (c),  the  solid  curves  are  the  raw  signal  strength  plus  noise, 
while  the  dashed  curves  are  noise  alone  at  a  slightly  detuned  frequency. 

(3)  In  (b)  and  (d),  the  final  estimates  of  the  TM  and  TE  field  strengths  are  ob¬ 
tained  using  Eqs.  (15)  and  (21). 

(4)  Profiles  are  drawn  only  where  both  the  signal -to -noise  and  the  recorded 
signal -to-calculated  cross-coupling  ratios  equal  or  exceed  6dB. 

In  part  (b)  of  Figures  10,  15,  and  16,  the  arrow  marks  the  signal  strength 
measured  on  the  ground  at  launch. 

A  significant  feature  seen  in  all  the  profiles  is  the  high  background  noise  over 
the  first  20km  of  altitude,  apparently  due  to  the  burning  solid  fuel  ARCAS  rocket 
motor.  This  noise  tends  to  mask  signals  in  that  altitude  range.  (The  S/N  ratio  of 
the  strongest  signal  recorded,  NAA  at  17.8kHz,  did  exceed  +6dB  on  the  TM  chan- 
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Figure  8.  Height -Gain  Profiles  for  16.4kHz  (JXZ):  (a)  Recorded  TM  Signal 
(Solid  Line)  and  Noise  (Dashed  Line);  (b)  TM  Signal  Compensated  for  Rocket 
Tilt;  (c)  Recorded  TE  Signal  (Solid  Line)  and  Noise  (Dashed  Line);  (d)  TE  Sig¬ 
nal  Compensated  for  Rocket  Tilt 


nel,  but  the  curves  for  the  other  transmitters  have  large  gaps  where  the  S/N  ratio 
was  too  small. )  Above  20  km,  where  the  rocket  motor  was  no  longer  firing,  the 
noise  is  20dB  or  so  smaller  on  both  channels,  and  many  signals  that  had  been  ob¬ 
scured  now  stand  out  clearly. 

i 

7.  DATA  ANALYSIS  AND  DISCUSSION 

1 

For  each  of  the  nine  transmitters.  Table  2  gives  solar  illumination,  ground 
path,  and  whether  TM  and/or  TE  signals  were  recorded  during  the  rocket  flight. 
The  first  eight  sources  listed  are  the  ground-based  TM  transmitters.  TM  signals 
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Figure  9.  Height -Gain  Profiles  for  17.  1  kHz  (UMS):  (a)  Recorded  TM  Signal 
(Solid  Line)  and  Noise  (Dashed  Line);  (b)  TM  Signal  Compensated  for  Rocket 
Tilt;  (c)  Recorded  TE  Signal  (Solid  Line)  and  Noise  (Dashed  Line);  (d)  TE  Sig¬ 
nal  Compensated  for  Rocket  Tilt 


were  received  from  each  of  these,  irrespective  of  propagation  conditions.  On  the 
other  hand,  converted  signals  (TE)  were  seen  from  only  four  of  these  sources, 
specifically  those  with  propagation  paths  that  were  entirely  under  nighttime  condi¬ 
tions  and  that  did  not  cross  over  the  Greenland  ice  cap.  For  those  transmitters, 
the  TE  fields  usually  were  somewhat  weaker  than  the  corresponding  TM  fields, 
although  at  aircraft  altitudes,  they  were  of  the  same  order  of  magnitude.  The 
other  four  sources,  for  which  no  converted  signals  were  measured,  had  propaga¬ 
tion  paths  involving  the  day /night  terminator  and  that  crossed  the  ice  cap  as  well. 

The  presence  of  TE  signals  from  ground-based  transmitters  (which  do  not  ra¬ 
diate  TE  polarization)  is  due  to  conversion  by  the  geomagnetic  field.  The  amount 
of  this  conversion  depends  in  part  on  the  effective  height  of  the  upper  (ionospheric) 
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Figure  11.  Height-Gain  Profiles  for  19.  05kHz  (GQD):  (a)  Recorded  TM  Signal 
(Solid  Line)  and  Noise  (Dashed  Line);  (b)  TM  Signal  Compensated  for  Rocket 
Tilt;  (c)  Recorded  TE  Signal  (Solid  Line)  and  Noise  (Dashed  Line);  (d)  TE  Sig¬ 
nal  Compensated  for  Rocket  Tilt 


larger  than  the  36  dB  >  ly  V/m  signal  observed  by  the  ground  receiver,  most  likely 
was  cross-coupled  TE  due  to  rocket  tilt;  thus,  no  TM  profile  is  shown  in  Figure 
16b.  This  absence  of  measured  TM  indicates  very  little  magneto-ionic  conversion 
of  TE  to  TM,  in  contrast  to  the  significant  TM  to  TE  conversion  observed  for 
ground-based  transmitters  with  similar  propagation  paths  (that  is,  nighttime,  no 
ice  cap. 

It  is  important  to  note  that  the  TE  signal  strengths  (converted  TM)  from  the 
ground-based  transmitters  in  the  nighttime  hemisphere  were  equal  to  or  were 
larger  than  the  TE  signals  received  from  TACAMO.  Although  consideration  must 
be  given  to  account  for  effects  due  to  differences  in  radiated  power,  frequency, 
range,  or  path  parameters  for  the  sources,  this  observation  has  strong  implica- 
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Figure  12.  Height -Gain  Profiles  for  19.  6  kHz  (GBZ):  (a)  Recorded  TM  Signal 
(Solid  Line)  and  Noise  (Dashed  Line);  (b)  TM  Signal  Compensated  for  Rocket 
Tilt;  (c)  Recorded  TE  Signal  (Solid  Line)  and  Noise  (Dashed  Line);  (d)  TE  Sig¬ 
nal  Compensated  for  Rocket  Tilt 


tions  in  the  assessment  of  coverage  of  airborne  systems  in  the  presence  of  other 
signals. 

A  feature  noticed  in  the  TM  height  profiles  is  the  tendency  for  those  which 

cross  the  ice  cap  (Figures  8,  9,  11,  12)  to  become  stronger  with  altitude,  while 

those  not  crossing  the  cap  (Figures  10,  13,  14,  15)  become  weaker  with  altitude. 

The  signals  from  the  two  transmitters  in  the  United  Kingdom  (Figures  11,  12)  show 

this  trend  quite  clearly,  the  amplitudes  increasing  by  more  than  lOdB  as  the  alti- 

fi 

tude  increases  from  20km  to  60km.  Waveguide  calculations  show  that  TM  profiles 


6.  Field,  E.C.,  Jr.,  and  Warber,  C.R.,  Pacific -Sierra  Research  Corp. ,  Los 
Angeles,  Calif.  ( 1985)  Private  communication. 
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Figure  13.  Height -Gain  Profiles  for  21.4kHz  (NSS):  (a)  Recorded  TM  Signal 
(Solid  Line)  and  Noise  (Dashed  Line);  (b)  TM  Signal  Compensated  for  Rocket 
Tilt;  (c)  Recorded  TE  Signal  (Solid  Line)  and  Noise  (Dashed  Line);  (d)  TE  Sig¬ 
nal  Compensated  for  Rocket  Tilt 


over  "Greenland  ice"  do  indeed  increase  with  altitude  (see  Figure  17).  It  appears 
that  the  dielectric-like  ice  cap,  in  effect,  distorts  the  waveguide  fields,  displacing 
the  TM  maximum  from  the  lower  boundary. 

The  rocket  was  launched  essentially  at  the  boundary  between  a  very  good  con¬ 
ductor  (sea  water)  and  a  very  poor  one  (the  polar  ice  cap).  The  transition  region 
is  very  abrupt,  the  conductivity  changing  many  orders  of  magnitude  within  a  few 
kilometers  (much  less  than  one  wavelength  at  20kHz)  of  the  launch  site.  The  data 
show  that  the  height  profiles  tend  to  retain  characteristics  corresponding  to  the 
ground  parameters  of  the  entire  propagation  path,  rather  than  immediately  trans¬ 
forming  to  new  profiles  as  they  cross  a  boundary.  Thus,  at  Thule  or  other  loca¬ 
tions  where  the  conductivity  changes  abruptly,  source  azimuth  is  more  important 
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Figure  14.  Height -Gain  Profiles  for  23,4kHz  (NPM):  (a)  Recorded  TM  Signal 
(Solid  Line)  and  Noise  (Dashed  Line);  (b)  TM  Signal  Compensated  for  Rocket 
Tilt;  (c)  Recorded  TE  Signal  (Solid  Line)  and  Noise  (Dashed  Line);  (d)  TE  Sig¬ 
nal  Compensated  for  Rocket  Tilt 


in  determining  signal  height -gain  than  is  the  ground  conductivity  directly  under  the 
receiver. 


8.  CONCLUSIONS 

This  rocket  experiment  has  provided  measured  height-gain  profiles  of  nine 
very  low  frequency  (VLF)  signals,  at  high  latitude  and  under  quiet  ionospheric  con 
ditions.  The  profiles  were  shown  to  depend  on  ground  conductivity  and  solar  illu¬ 
mination  conditions  on  the  propagation  paths.  Conversion  of  TM  to  TE  was  quite 
evident  for  paths  under  total  nighttime  conditions,  while  no  conversion  was  ob- 
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Figure  16.  Height -Gain  Profiles  for  26.  1kHz  (TACAMO):  (a)  Recorded  TM  Signal 
(Solid  Line)  and  Noise  (Dashed  Line);  (b)  TM  Signal  Compensated  for  Rocket 
Tilt;  (c)  Recorded  TE  Signal  (Solid  Line)  and  Noise  (Dashed  Line);  (d)  TE  Sig¬ 
nal  Compensated  for  Rocket  Tilt 
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Frequency 

Source 

Path 

Illumination 

Path  Over 
Ice  Cap 

Received 
TM  TE 

16.4 

Aldra,  Norway 

day /night 

yes 

yes 

no 

17.  1 

Gorki,  U.  S.  S.  R. 

day /night 

yes 

yes 

no 

17.8 

Cutler,  Maine 

night 

no 

yes 

yes 

19.05 

U.K. 

day /night 

yes 

yes 

no 

19.6 

Criggion,  U.K. 

day  /  night 

yes 

yes 

no 

21.4 

Annapolis,  Md. 

night 

no 

yes 

yes 

23.4 

Haiku,  Hawaii 

night 

no 

yes 

yes 

24.8 

Jim  Creek,  Wash. 

night 

no 

yes 

yes 

26.  1 

TACAMO 

night 

no 

no 

yes 

Table  2.  Propagation  Conditions 
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Figure  17.  Calculated  Height-Gain  Profiles  for  GBZ  to  Thule 
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